Numerical solution of the Optimal Transportation problem using the Monge-Ampère equation
نویسندگان
چکیده
A numerical method for the solution of the elliptic MongeAmpère Partial Differential Equation, with boundary conditions corresponding to the Optimal Transportation (OT) problem is presented. A local representation of the OT boundary conditions is combined with a finite difference scheme for the Monge-Ampère equation. Newton’s method is implemented leading to a fast solver, comparable to solving the Laplace equation on the same grid several times. Theoretical justification for the method is given by a convergence proof in the companion paper [BFO12]. In this paper, the algorithm is modified to a simpler compact stencil implementation and details of the implementation are given. Solutions are computed with densities supported on non-convex and disconnected domains. Computational examples demonstrate robust performance on singular solutions and fast computational times.
منابع مشابه
Numerical solution of the second boundary value problem for the Elliptic Monge-Ampère equation
This paper introduces a numerical method for the solution of the nonlinear elliptic Monge-Ampère equation. The boundary conditions correspond to the optimal transportation of measures supported on two domains, where one of these sets is convex. The new challenge is implementing the boundary conditions, which are implicit and non-local. These boundary conditions are reformulated as a nonlinear H...
متن کاملA numerical method for the elliptic Monge-Ampère equation with transport boundary conditions
The problem of optimal mass transport arises in numerous applications including image registration, mesh generation, reflector design, and astrophysics. One approach to solving this problem is via the Monge-Ampère equation. While recent years have seen much work in the development of numerical methods for solving this equation, very little has been done on the implementation of the transport bo...
متن کاملA Least-Squares Method for Optimal Transport Using the Monge-Ampère Equation
In this article we introduce a novel numerical method to solve the problem of optimal mass transport and the related elliptic Monge-Ampère equation. It is one of the few numerical algorithms capable of solving this problem efficiently with the proper boundary conditions. It scales well with the grid size and has the additional advantage that the target domain may be non-convex. We present the m...
متن کاملSolving the Monge-Ampère Equations for the Inverse Reflector Problem
The inverse reflector problem arises in geometrical nonimaging optics: Given a light source and a target, the question is how to design a reflecting free-form surface such that a desired light density distribution is generated on the target, e.g., a projected image on a screen. This optical problem can mathematically be understood as a problem of optimal transport and equivalently be expressed ...
متن کاملMesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation
Article history: Received 12 June 2015 Received in revised form 7 December 2015 Accepted 9 December 2015 Available online 15 December 2015
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 260 شماره
صفحات -
تاریخ انتشار 2014